YIPENG LIU: ROBUST CAPON BEAMFORMING VIA SHAPING BEAM PATTHiN 



Robust Capon Beamforming via Shaping Beam 

Pattern 

Yipeng Liu 



Abstract 

High sidelobe level and direction of arrival (DOA) estimation sensitivity are two major disadvantages of the 
Capon beamforming. To deal with these problems, this paper gives an overview of a series of robust Capon 
beamforming methods via shaping beam pattern, including sparse Capon beamforming, weighted sparse Capon 
i beamforming, mixed norm based Capon beamforming, total variation minimization based Capon beamforming, 

mainlobe-to-sidelobe power ratio maximization based Capon beamforming. With these additional structure-inducing 
H ' constraints, the sidelobe is suppressed, and the robustness against DOA mismatch is improved too. Simulations 

' show that the obtained beamformers outperform the standard Capon beamformer. 
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Index Terms 

array signal processing, robust Capon beamforming, beam pattern shaping, sidelobe suppression, direction of 
arrival (DOA) mismatch, sparse constraint, cosparse constraint. 



c/3 ' 1. Introduction 

o ■ 

' — Beamforming is a conventional technique for signal's emission and reception in some certain directions 
in sensor array system [HJ [)2] [|3J. To selectively receive and transmit in space, beam pattern can be 
>■ synthesized either adaptively or deterministically. Comparing with the omnidirectional sensor system, 
beamforming system exploits the signal's spatial dimension to enhance the signal quality. It has a wide 

^ range of applications in the field of radar, sonar, acoustics, astronomy, seismology, communications, and 

^ medical imaging. 

CN The non-adaptive (data-independent) beamforming consists of the delay-and-sum method and a variety 
^ of weighting based sidelobe control methods; while adaptive (data-dependent) beamforming solves a 
,—1 optimization problem with a data-driven performance function to get a set of array weighting vector. The 
j>! data-independent beamformer uses a set of pre-defined weights to linearly combine the transmitted/received 
signals in different sensors, which only uses the information about spatial position of the sensor and signal- 
of-interest (SOI). Generally data-dependent beamforming will further exploit some characteristics of the 
transmitted/received signal, to suppress the interference and noise in the non-interested directions. 

The Capon beamformer is one of the most popular adaptive beamforming systems. It minimizes the 
array output power while subjecting to the linear constraint that the SOI does not suffer from any distortion 
by adaptive selection of the weighting vector. The Capon beamformer has better resolution and much better 
interference rejection capability than the data-independent beamformer. However, its high sidelobe level 
and the SOI steering vector uncertainty due to differences between the assumed signal arrival angle and 
the true arrival angle would seriously degenerate the performance in the presence of environment noise 
and interferences [j4j [|5l. 

With a spherical uncertainty set being introduced, doubly constrained robust (DCR) Capon beamformer 
is obtained with the increased robustness against direction of arrival (DOA) mismatch flU. The diagonal 
loading method is one of the most popular ways to deal with the uncertainty in steering vector in robust 
Capon beamforming (RCB) [HI [|71. But its main drawback is that the diagonal loading factor is not 
convenient to determine. Eigen-space-based beamforming approach [8], which can only be used to the 
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point signal source and high signal-to-noise-ratio (SNR) cases, uses the projection of the presumed signal 
steering vector onto the sample signal-plus-interference subspace instead of the presumed signal steering 
vector. Many works have been done to enhance the robustness against the SOI steering vector uncertainty, 
to accelerate the convergence, etc [1]. 

In this paper, different from the previous RCB, a class of new methods to obtain the robustness by 
shaping the beam pattern are summarized. In [9|, sparse constraint is used to encourage sparse distribution 
of the array gains. In [10], weighting is used to make the distribution sparser. They are the first two papers 
that use the beam pattern shaping constraints. But the structure information of the beam pattern is not 
fully exploited. In order to better match to the ideal beam pattem, several other constraints are proposed, 
including the mixed norm based constraint [[TTI . total variation minimization (TVM) constraint [fT2l [[T3l . 
and the mainlobe to sidelobe power ratio maximization constraint [|14il . These beam pattem shaping 
constraints use different ways to encourage dense distribution of the array gains in the mainlobe and 
sparse distribution of the array gain in the sidelobe. i.e., nearly all the large gains are accumulated in the 
vicinity of the mainlobe, while the smaller ones are in the sidelobes. Thereby, when there is uncertainty 
in the estimated direction of SOI, the array gain in the real direction of SOI, which is often near the 
estimated one, can be large enough to keep out large SOI power loss. At the same time, the interferences 
and noise is largely blocked, because nearly all the array gains in the sidelobes are much smaller. Numerical 
experiments show that the beam pattem shaping constraint based robust Capon beamforming outperforms 
standard one. Different kinds of beam pattem shaping constraints have their own advantages. 

II. Signal Model 

As it shows in Fig. [U the signal impinging into a uniform linear array (ULA) with M antennas can be 
represented by an M-by-1 vector [HI: 

J 

x(k) = simiOo) + J]/3j(k)a(ej) + n(k) (1) 

j=i 

where k is the index of snapshot, J is the number of interferences, s{k) and /3j{k), j = 1, ... , J, are the 
amplitudes of the SOI and interferences at k, respectively, dj (for / = 0, 1, ... , 7) are the values of DOA 
of the SOI and interferences. n(k) is the additive white Gaussian noise (AWGN) vector at time instant k. 
a(6j), j = 1, ... , J is the array steering vector, whose m-th element is 

[aiOj)]^ = exp |/(m - 1)^ sin 6^ (2) 

where / = V-T, d is the distance between two adjacent sensors and A is the wavelength of the SOI. The 
output signal of a beamformer at the time instant k can be formulated as: 

J 

y{k) = w"xik) = s{k)w"aieo) + _^/3//t)w^a(0y) + w^n(yt) (3) 

where w is the M-by-1 complex- valued weighting coefficients of beamformer. 

III. Classical Capon beamforming 

In classical adaptive beamformer, it enforces array gain in the estimated direction to be constant 
w^a(0o) ~ 1 and minimize the array output power of interference and noise w^n(^) to make the array 
output is approximately equal to the SOI. Here the array response gain in the direction is w^a(6^o), and 
the noise and interference gains in the array output is w^R„w, wherein R„ = E^nn^j, and E denotes the 
expectation operator. Assuming a(6o) and R„ are known in advance, the optimization model to produce 
the optimal array weighting vector is: 

min(w^R„w), s. t. w^a(0()) = 1 (4) 
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where R„ and Oq axe assumed to be known exactly in advance. However, it is not always in this case. 
The Capon beamformer is proposed to replace the noise and interference power R„ with the covariance 
of the received signals R^. The estimated covariance of the received signals R^ e qMxm obtained 
by processing multiple snapshots received by multiple sensors, i.e.: 

1 ^ 

R. = T7 X '^OVO') (5) 

j=k-M+l 

Thus, the Capon beamforming can be formulated as: 

Wmvdr = arg min (w^Rjcw) , s. t. w^a{ao) = 1 (6) 

w 

where w^R^w is the minimum covariance constraint, w^a(ao) = 1 is the distortionless constraint of the 
SOI. That why the Capon beamforming is called minimum variance distortionless response (MVDR) 
beamforming too; ao is the estimated direction of SOI, which cannot be exactly equal to the real direction 
of SOI ^0 because of estimation error. Solving the optimization model of MVDR beamforming can give 
a close form solution for the optimal weighting vector as: 

R;'afa)) 

y^MVDR - —777 , I / ) 

The difference between (H) and Q is that the MVDR beamforming minimizes the power of SOI, 
interference and noise, and meanwhile its distortionless constraint can keep there is little loss in the 
direction of SOI, which guarantees the array gain in the direction of SOI is not degenerated. 

IV. Capon beamforming with shaped beam pattern 

A. Sparse Capon beamforming 

In [9], a sparse constraint is incorporated into the MVDR beamformer to enforce sparse distribution in 
the whole beam pattern. It can reduce the number of nonzero array gains and suppress the sidelobe level 
of the classical MVDR beamformer. The sparse constraint (SC) based improved MVDR beamformer can 
be formulated as: 

wsc = arg min (w^R.vW + 71 ||w^A||^) , s. t. w^a(0o) = 1 (8) 

where the non-negative scalar 71 is the parameter that makes balance of the minimum variance constraint 
and the sparse constraint, ||x||p = (2r l-^iD^^'' is the Lp norm of a vector x, the M-hy-N A is the array 
manifold with a„s ( n = I, 2, ... , N) being the sampled angles ranging from -90° to 90°, which covers 
all the steering vectors for signals impinging from all possible angles, i.e. 



A = 



1 ••• 1 

expOVi) ■•■ expOVw) 

_ exp (j (M - 1) ^1) • • • exp (j (M -1)(Pn) 



(9) 



2nd 

(p„ = —-sman, for n= 1,2,..., N (10) 
A 

Different from previous sparse constraint in synthesis form, i.e. the sparsity is obtained by decomposing 
the estimated variables into a synthesis dictionary and a sparse vector, the used sparsity constraint here is 
in analysis, i.e. the enforced sparsity in the Lp norm constraint is obtained by multiplying the estimated 
variables with a synthesis dictionary. Recently this kind of sparsity in analysis form is called cosparsity 
too Ml. 

Considering that when < p < 1, the minimization of the Lp norm is not convex, = 1 is chosen to 
make the optimization model for sparse beamformer be convex, i.e.: 
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Wsc = arg min (w^R;tW + yi ||w^A|U , s. t. w^a(Q;o) =1 (11) 

w 

The sparse Capon beamforming (fTTT l is a convex programming, and can be solved efficiently [[T5l . 

B. Weighted sparse Capon beamforming 

We consider the M-by-l vector, x{k), as a snapshot of the received signal at time instant k. If we collect 
the snapshots of ^ > 1) different time instants in a matrix, then we can have an M-by-^ data matrix 

as 

X= [x(l),x(2),...,x(i^)] (12) 

It has been shown that the cross-correlation of the steering matrix A and the received data matrix X 
coarsely represents the a posteriori spatial distribution of interfering signals [17J. We use this property to 
define a weighted sparse constraint for further suppressing sidelobe level of the beam pattem. As a result, 
the weighting vector of the beamformer with a weighted sparse constraint is given by 

Wivsc = arg min (w^R;cW + 72 ||w^AQ||''j , s. t. w^a(0o) = 1 (13) 

where the N-hy-N matrix Q = diag[S N M (A^X^] serves as a weighting matrix, and SNM(a"X^ is an 
A^-by-1 vector containing elements the squared normalized mean value of each row of the N-hy-K matrix 
A^X ifTOl [fTTl . and 72 is the weighting factor balancing the minimum variance constraint and the beam 
pattem shaping constraint. Comparing ([131 ) with (fTTT i. we can see that the matrix Q in (fT3l) provides 
additional weighting on the sparse constraint, in accordance with the DOA distribution of interfering 
signals. More specifically, the larger the probability of interference arriving in a certain direction, the 
larger the weight applied on the sparse constraint in the corresponding direction. 

The optimal weighting vector indicated by (fT3T ) can be found by using an adaptive iteration algorithm 
[fTOl . When p = I, a series of algorithms for convex programming, can be used to solve (fT3l) efficiently 
[fTSl . We also observe that (fTT)) can be considered as a special case of (fT3l) . in terms that (fT3l) reduces to 
(fTT|) when Q = I, corresponding to the case of equal weighting in every direction. 

C. Mixed norm based Capon beamforming 

In (fTTj) . the added constraint encourages sparse distribution in all the possible values of DOA from 
-90° to 90°, no matter whether the array gains are in the mainlobe or the sidelobe. However, the array 
gains are not in this kind of conventional sparse distribution. The array gains are densely distributed in 
the mainlobe and sparsely distributed in the sidelobe. To exploit this more detailed structural sparsity 
information to improve the performance, a mixed norm constraint with two kinds of norms on mainlobe 
and sidelobe respectively can be added to the Capon beamformer. The new beamformer can be formulated 
as [fTTl: 

Wmnc = argmin [w^R.vW + 73 (||w^Am|L + llw^A^IIJ] , s. t. w^a(6io) = 1 (14) 

w 

where 73 is the weighting parameter. Am and A5 contains the vectors from the steering matrix A, A^ is 
composed of + 1 steering vectors corresponding to the mainlobe; and A5 is constituted with the steering 
vectors in A corresponding to the sidelobe. The product w^A^ indicates array gains of the mainlobe in 
the beam pattem, and w^As indicates array gains of the sidelobe. b is an integer in represention of the 
bound between the mainlobe and the sidelobe. T minimization of Hw^A^H^ enforces the sparse distribution 
of the array gains in the sidelobe, and the minimization of ||w^Am||^ results in dense distribution of the 
array gains in the mainlobe. 

To illustrate why different kinds of norms encourage different kinds of distributions, a simple 2- 
dimensional geometry illustration is given in Fig. |2] The minimization of the LI norm, the L2 norm and 
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the Loo norm of a two dimensional vector can be represented as a smallest rhombus, a smallest circular 
and a smallest exact square, respectively. The distortionless constraint is a fixed line in the plane. The 
optimal solution for the minimization of the norm ||w^A||^, (p = 1, 2, oo) subjecting to the distortionless 
constraint w^a(0()) = 1 would be the tangent point of the line and the curve (rhombus, circular, square). 
Therefore, the LI norm minimization results in two entries of the solution which have quite different 
absolute values with high probability; and the Loo norm minimization gives birth to two entries of the 
solution which have considerable similar absolute values with high probability. 

With this mixed norm constraint added, the obtained beam pattern would be sharped like this: most 
of the significant array gains are located in the mainlobe area; and the rest trivial array gains are in the 
sidelobe area. Since most of entries in the mainlobe are non-trivial, the array gains in the mismatched 
angle are also significant, which avoids seriously degenerate the performance. In addition. Since most of 
entries in the sidelobe are trivial, nearly all of the array gains in the angles of interferences and noise are 
very small, which reduce the power of interferences and background noise in the output signal. 

Here we call (fT4l) the mixed norm (MN) based Capon (MNC) beamformer, and its solution can be 
given by some efficient algorithms IfTSl . 



D. Total variation based Capon beamforming 

The sparse constraint encourages sparse distribution for all the array gains w^A for all the possible 
values of DOA from -90° to 90°. However, the array gains in the mainlobe are not sparse, but dense as a 
solid block. To improve the performance with a more suitable constraint on the beam pattern, the LI norm 
minimization based sparse constraint is only added on the sidelobe, and a TVM restricts for the entire 
beam pattern [|T2]| [fT3l . The TVM and sparse constraint based Capon beamformer can be formulated as 

m-. 

Wtvms = arg min |w^R;cW + 74 |x; ||d,(w^A)^||^ + Hw^A^] 



s. t. w^a(Qro) = 1 



(15) 



where 



D, 



D,5 



(16) 



Di f and D,b are the i-th order forward and backward differential matrices. / is the total number 
of differential matrix D,; As is constituted with the sidelobe steering vectors in A. The product w^As 
indicates array gains of the sidelobe. 74 is the weighting factor controlling the TVM constraint and the 
sparse constraint. Since the objective function of the proposed beamformer (fTST ) is convex; the optimal 
Wtvms can also be solved out by convex programming software [flSl . 

In the beam pattern shaping constraint 2 D,(w^A) + llw^Asll the first term is the TVM which 

discourages the large fluctuation in the beam pattern. It results in the high array gains accumulated in 
the mainlobe and the small trivial ones gathered in the sidelobe. In addition, the sparse beam pattern 
constraint is modified to be the second term to further suppress the sidelobe level. As the new constraint 
in (fT5l) fits the desired beam pattern better, the performance would be enhanced. 



E. Mainlobe-to-sidelobe power ratio maximization based Capon beamforming 

In the perspective of the beam pattern, it is observed from the Capon beamformer that there is only an 
explicit constraint on the desired DOA, i.e. w^a(0o) = 1. while no constraint is put onto the interference and 
background noise. To repair this drawback, we propose the following cost function with a regularization 
term, which forces maximization of the MSPR [[T4l : 
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WM5PR = arg mm < w^R.^w + js" 



s. t. w^a(0o) = 1 



^"^m\\1 I (17) 



where 75 is the weighting factor balancing the minimum variance constraint and the MSPR maximization 
constraint. 

Physically w"a{6) is the array gain in the signal direction 9. The product w^Aj indicates array 
gains of the sidelobe; and the product w^A^ indicates array gains of the mainlobe. With the term 

IIw^AsII^^IIw^AmIIj minimized, the MSPR based Capon beamformer (fTSl) can perform sidelobe min- 
imization and mainlobe maximization simultaneously. But unfortunately it is not convex. To make it 
convex, we relax the MSPR constraint and obtain a new beamformer as 



y^RMSPR = arg min jw^R.vW + 75 (||w^Am||2 " l) + Hw^A^ || 

s. t. w"a(eo) = 1 



(18) 



The splitting of the matrix A into Am and Am helps. The newly added relaxed MSPR (RMSPR) 
term ^||w^Am||2 - ij + ||w^A5||2 is convex. It is minimized to minimize the sidelobe power Hw^A^H^ 

and the approximation error ^||w^Am||2 ~ ^) • '^^^ ^i'^^ °f '•^ reshaping the Capon beam 

pattern with constraint on both mainlobe and sidelobe. The mainlobe power and sidelobe power are 

II 1 12 

constrained separately in the optimization model (fTST ). It is no doubt that the minimization of Hw^AsH^ 
gives smaller power of sidelobe. The minimization ^Hw^A^H^ - ij makes the power of mainlobe be a 

constant approximately, i.e. for ^Hw^A^H^ - ij < ^, if ^ is very small, the power of mainlobe ||w^Am||2 
can be approximately constant. Specially if ^ = , the mainlobe power is strictly constant. Combining these 
two constraints for beam pattern shaping, it forces the power of mainlobe Hw^A^H^ to be an approximately 

constant while making the power in sidelobe Hw^A^H^ as small as possible. With the power in mainlobe 
being a constant and the power in sidelobe being minimized, the MSPR can be maximized in a relaxed 

II 112/11 ii2 

way. Thus, the relaxed form of the maximization of MSPR Hw^A^ly ||w^Am||2 would be achieved. 

V. Numerical experiments 

Numerical experiments are employed to demonstrate the performance of different beam pattern shaping 
constraint based Capon beamforming methods and the standard Capon beamforming. a ULA with 8 half- 
wavelength spaced antennas is considered. The AWGN at each sensor is assumed spatially uncorrelated. 
The DOA of the SOI is set to be 0°, and the DOAs of three interfering signals are set to be -30°, 30°, 
and 70°, respectively. The SNR is set to be 10 dB, and the interference to noise ratios (INRs) are assumed 
to be 20 dB, 20 dB, and 40 dB in -30°, 30°, and 70°, respectively. 100 snapshots are used for each 
simulation. The matrix A consists of all steering vectors ranging in [-90°, 90°] with the sampling interval 
of 1°. Without loss of generality, p is set to be 1; b is set to be 15; and y,, /=1,2,...,5 are all set to get 
the optimal performance. 

To quantify the performance evaluation of different beamformers, the signal to interference noise ratio 
(SINR) of beamformer's output can be defined as 

cr>^a(0o)a^(eo)w 
SINR = 7y " ^ " ^ (19) 



w^j^Z cr2a(e,)a^(e,) + Q)w 



YIPENG LIU; ROBUST CAPON BEAMFORMING VIA SHAPING BEAM PATTHiN 



7 



where cr^ and cr^ are the variances of the SOI and the j-th interference, and Q is a diagonal matrix whose 
diagonal elements are the variances of the noise. 

When the estimated DOA of the SOI and the real one are exactly the same, i.e. Oq = ao = 0°, in 
1000 times Monte Carlo simulations. Fig. [3] gives the beam pattems of the Capon beamforming, sparse 
beamforming and weighting sparse Capon beamforming; Fig. |4] gives the beam pattems of the Capon 
beamforming, the sparse Capon beamforming and the mixed norm based Capon beamforming; Fig. |5]gives 
the beam pattems of the Capon beamforming, the sparse Capon beamforming and the TVM based Capon 
beamforming. Fig. |6] gives the beam pattems of the Capon beamforming, the sparse Capon beamforming 
and the MSPR based Capon beamforming. Each beam pattem is normalized with its L2 norm, i.e. the 
power of the array gains is 1. Obviously it can be seen from the figures that the RCBs with pattem 
shaping constraints outperform the standard one. Compared with the sparse Capon beamforming, the 
weighted sparse Capon beamforming, the mixed norm based Capon beamforming, the TVM based Capon 
beamforming and the MSPR based Capon beamforming have much lower array gains in the directions 
of interferences (-30°, 30°, and 70°). They have better nulling performance for interference suppression, 
and the SINR of the array output signals. In 1000 times Monte Carlo simulations, the average SINR of 
the beamformers are 2.3027 dB, 4.3178 dB, 5.0162 dB, 5.8119 dB, 5.8563 dB and 6.5224 dB. 

In practice, the high sensitivity of angle mismatch is one of the main disadvantages of Capon beam- 
forming. In standard Capon beamforming, when the estimated DOA of SOI is quite different from the 
real one, i.e. 6o 4^ ao, the distortionless constraint of the Capon beamforming is for the direction, which 
allows the signal from the direction ao is fully received. However, the actual signal is from the direction 
aoj the other constraint of the Capon beamforming, i.e. the minimization of the output power of the 
array, would regard the SOI from direction ao as an interference. The minimum variance distortionless 
constraints would make a nulling in the direction of SOI, which would greatly decrease the SINR of the 
output signals. When there is a 3° angle mismatch (^o = 3° and ao = 0°), the beam pattems are shown in 
Fig. |7] and Fig. [lOl Suppressing the sidelobe level and obtaining deeper nuUings to avoid interferences. 
Fig. |71 Fig. [8l Fig. |9] and Fig. [10] show that beam pattem shaping constraint based Capon beamforming 
has better robustness than standard Capon beamforming. Obviously we can see that there is a nulling in 
the direction of SOI. However, the array gains in the real DOA of SOI have almost the same level of the 
estimation DOA. Similarly, compared sparse Capon beamforming, weighted sparse Capon beamforming, 
mixed norm based Capon beamforming, and TVM based Capon beamforming, and MSPR based Capon 
beamforming have lower nulling for avoiding interference. When there is angle mismatch, as denoted in 
table m in 1000 times' Monte Carlo simulations, the values of SINR of each beamformers' output signal 
are: 0.0003 dB, 3.1903 dB, 3.8013 dB, 4.6836 dB, 4.8667 dB and 3.8402 dB. 

In the above discussion, the SINR values show the performance of the sparse Capon beamforming is the 
worst of all the beam pattem shaping constraints based Capon beamforming, but it has some advantages 
in some other applications. Comparing the beam pattems of sparse Capon beamforming and other robust 
Capon beamforming methods, the beam pattem of sparse Capon beamforming has a narrow mainlobe, 
which is preferred in radar applications. 

In addition, in the presence of 3 degrees mismatch, the SINR of MSPR constraint based Capon 
beamforming is less than the values of SINR of mixed norm based Capon beamforming and TVM 
based Capon beamforming. However, in the case of no angle mismatch, MSPR constraint based Capon 
beamforming has the best SINR performance. Thus we can see that MSPR constraint based Capon 
beamforming has best performance of suppressing sidelobe, but its robustness against the angle mismatch 
is inferior to the mixed norm based Capon beamforming and TVM based Capon beamforming. 

To sum up, a large performance improvement can be obtained by adding the proposed beam pattem 
shaping constraints with respect to standard Capon beamforming. Meanwhile, Comparing the beamforming 
methods in this paper, each methods have its own characteristics. In certain environments, each has some 
unique performance advantages. 
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VI. Conclusion 

This paper begins with a brief introduction of the Capon beamforming research and existing problems. 

Then, the problems of the standard Capon beamforming are given: the high sensitivity of the DOA 
mismatch and the high sidelobe level. A set of novel methods are summarized to deal with these two 
issues. The methods used here improve the beamforming performance by adding the array gain distribution 
encouragement constraints. In order to get an ideal beam pattern distribution. This paper has discussed 
the sparse constraint, the weighted sparse constraint, the mixed norm constraint, TVM constraint, and the 
MSPR maximization constraint. Numerical experiments show that the performance of Capon beamforming 
is significantly improved by using beam shaping constraints. It preferably overcomes the high sensitivity of 
angle mismatch problem and high sidelobe level. In addition, the performance of different beam-shaping 
constraints are compared. 

In the future, the beam pattern shaping constraint can be combined with other robust Capon beamforming 
techniques, such as ellipsoid methods, diagonal loading methods, moreover, it is natural to generalize the 
beam pattern shaping constraints to the 3-D beamforming. 
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TABLE I 

The average values of SINR of the beamformers with and without 3° DOA mismatch (dB). 





Capon 


sparse Capon 


weighted sparse Capon 


mixed norm 


TVM 


MSPR 


no DOA mismatch 


0.0281 


0.1492 


0.0527 


0.2146 


0.1012 


0.2047 


3° DOA mismatch 


0.0231 


0.1283 


0.0577 


0.1642 


0.1562 


0.1860 




Fig. 1. Array structure of the beamformer. 
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Fig. 2. The geometric illustration of the Lp noiin minimization when p = 0, 1,2, oo. 
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Fig. 3. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and weighted sparse beamforming, without any 
mismatch between the estimated DOA of SOI and the real one. 
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Fig. 4. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and mixed norm shaped Capon beamforming, 
without any mismatch between the estimated DOA of SOI and the real one. 
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Fig. 5. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and TVM shaped Capon beamforming, without 
any mismatch between the estimated DOA of SOI and the real one. 
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Fig. 6. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and MSPR shaped Capon beamforming, without 
any mismatch between the estimated DOA of SOI and the real one. 
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Fig. 7. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and weighted sparse Capon beamforming, with 
3° mismatch between the estimated DOA of SOI and the real one. 
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- The MN based Capon beamformer 
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Fig. 8. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and mixed norm shaped Capon beamforming, with 
3° mismatch between the estimated DOA of SOI and the real one. 
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Fig. 9. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and TVM shaped Capon beamforming, with 3° 
mismatch between the estimated DOA of SOI and the real one. 
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Fig. 10. The beam patterns of standard Capon beamforming, sparse Capon beamforming, and MSPR shaped Capon beamforming, with 
mismatch between the estimated DOA of SOI and the real one. 



